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Abstract. We describe an automated method for assigning the most likely physical 
parameters to the components of an eclipsing binary (EB) , using only its photometric 
light curve and combined color. In traditional methods (e.g. WD and EBOP) one 
attempts to optimize a multi-parameter model over many iterations, so as to min- 
imize the chi-squared value. We suggest an alternative method, where one selects 
pairs of coeval stars from a set of theoretical stellar models, and compares their 
simulated light curves and combined colors with the observations. This approach 
greatly reduces the EB parameter-space over which one needs to search, and allows 
one to determine the components' masses, radii and absolute magnitudes, without 
spectroscopic data. We have implemented this method in an automated program us- 
ing published theoretical isochrones and limb-darkening coefficients. Since it is easy 
to automate, this method lends itself to systematic analyses of datasets consisting of 
photometric time series of large numbers of stars, such as those produced by OGLE, 
MACHO, TrES, HAT, and many others surveys. 



1. Introduction 

Eclipsing double-lined spectroscopic binaries provide the only method 
by which both the masses and radii of stars can be estimated with- 
out having to resolve spatially the binary or rely on astrophysical 
assumptions. Despite the large variety of models and parameter-fitting 
implementations (e.g. WD and EBOP), their underlying methodology 
is essentially the same. Photometric data provides the light curve of 
the EB, and spectroscopic data provide the radial velocities of its 
components. The depth and shape of the light curve eclipses constrain 
the components' brightness and fractional radii, while the radial ve- 
locity sets the length scale of the system. In order to characterize 
fully the components of the binary, one needs to combine all of this 
information. Unfortunately, only a small fraction of all binaries eclipse, 
and spectroscopy with sufficient resolution can be performed only for 
bright stars. The intersection of these two groups leaves a pitifully small 
number of stars. 

In the past decade, there has been a dramatic growth in the number 
of stars with high-quality, multi-epoch, photometric data. This has been 
due to major advances in both CCD detectors and the implementation 
of image-difference analysis techniques (Crotts, 1992; Alard et al., 1998; 
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Alard, 2000), which enables simultaneous photometric measurements 
of tens of thousands of stars in a single exposure. Today, there are 
many millions of light curves available from a variety of surveys, such 
as OGLE, Udalski et al. (1994); MACHO, Alcock et al. (1998); TrES, 
Alonso et al. (2004); and HAT, Bakos et al. (2004). But there has not 
been a corresponding growth in the quantity of spectroscopic data, nor 
is this likely to occur in the near future. Thus, the number of fully- 
characterized EBs has not grown significantly. In recent years there 
has been a growing effort to mine the wealth of available photometric 
data, by employing simplified EE models in the absence of spectroscopic 
observations (Wyithe et al., 2001, 2002; Devor, 2004, 2005). 

In this paper we present a novel approach, which utilizes theoretical 
models of stellar properties to estimate the orbital parameters as well as 
the masses, radii, and absolute magnitudes of the stars, while requiring 
only a photometric light curve and an estimate of the binary's combined 
color. This approach can be used to characterize quickly large numbers 
of eclipsing binaries, however it is not sufRcicnt to improve stellar mod- 
els since underlying isochrones must be assumed. We have created two 
implementations of this idea. The first, which we have named MECI- 
express, and is described in section 2, is a "quick and dirty" program 
that is designed as a simple extension to the Detached Eclipsing Binary 
Light curve (DEBiL) fitter (Devor, 2004, 2005). The second, which we 
have named MECI, and is described in section 3, is considerably more 
accurate, but also more computationally demanding. The source code 
for both MECI-express and MECI will be provided upon request. 



2. Express Method for Eclipsing Component Identification 

(MECI-express) 

The primary application of MECI-express is to identify the stellar 
components of a given EB. It operates after a conventional EB model- 
fitting program has already analyzed the given EB's light curve. In 
our implementation, we chose to employ DEBiL (Devor, 2004, 2005) 
since it is simple, fast, and fully automated. The fitted parameters 
are the orbital period (P), the apparent magnitudes {magi^2), and the 
fractional radii (ri^2) of the binary components. A fractional radius is 
defined as the radius (i?i,2) divided by the sum of the components' 
semimajor axes (a). In MECI-express we iterate through a large group 
of MK spectral type pairings, to each of which we associate typical 
stellar parameters (Cox, 2000). These stellar parameters are the masses 
(-^1,2), the radii (i?i,2), and the absolute magnitudes (Magi^i) of the 
binary components. If the assumed values of the stellar parameters 
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match the true values, then the stellar and fitted parameters should 
obey to the following equations: 



G{Mi + M2) 



P\t (1) 



G(Mi + M2) - ^ 
Magi — Mag2 = magi — mag2 (3) 

We also may have additional constraints from the observed out- 
of-eclipse combined colors of the system. For example, in the case of 
OGLE II targets, we have the estimated V-I color: 



Magv — Magi = magy — magi (4) 

We assume that the color has been corrected for reddening and 
that no systematic errors are present, so any inequalities would be due 
to an incorrect choice for the component pairing. The likelihood of 
each pairing is assessed by calculating the difference between the left- 
hand-side (stellar parameters) and right-hand-side (fitted parameters) 
of each equation. These differences are divided by their uncertainties, 
and added in quadrature. The pairing with the smallest sum is deemed 
the most likely pairing. For each given EB light curve, MECI-express 
returns the list of the top ranked (most likely) binary pairings, with 
their corresponding sums. MECI-express can also be used to create a 
contour plot of the probability distribution for all pairings. We illustrate 
an example of individual MECI-express components in Figures l.a-c, 
which are then combined to create the result shown in Figure 2. a. 




05 BO B5 AO A5 FO F5 GO G5 KO K5 WO M2 M5 05 BO B5 AO A5 FO F5 GO G5 KO K5 MO U2 M5 05 BO B5 AO A5 FO F5 GO G5 KO K5 MO M2 W5 

Storl type Storl type Storl type 



a. Constraints from (cq. 1) b. Constraints from (cq. 2) c. Constraints from (eq. 3) 

Figure 1, Contour plots of the absolute difference between the left-hand-side and the 
right-hand-side of each equation, divided by its uncertainty, as applied to the WW 
Camelopardalis light curve (Lacy et al., 2002). Adding these results in quadrature, 
produces the likelihood plot shown in Figure 2. a. 
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a. MECI-express likelihood plot 



b. MECI likelihood plot (age = 0.6 Gyr) 



Figure 2. A comparison of the MECI-express (left) and MECI (right) likelihood 
contour plots for WW Carnolopardalis. The value of the contours arc described in 
the body of the text. The Eisterisk marks the solution of Lacy et al. (2002). 



3. Method for Eclipsing Component Identification (MECI) 

MECI was developed to improve significantly upon the accuracy of 
MECI-express (see Table 1). This was done as follows: We replaced the 
use of spectral types with the more fundamental (and continuous) quan- 
tities of mass and age. Furthermore, in MECI we assume that the two 
binary components are coeval, thus replacing the 2-dimensional spectral 
type - spectral type grid, with a 3-dimensional mass-mass-age grid. 
Finally, we no longer rely on parameter fits of the components' apparent 
magnitudes and fractional radii directly from the light curve, which are 
often very uncertain, nor do we assume constant limb-darkening coeffi- 
cients. Instead, we interpolate these values for the given mass-mass-age 
pairing, from precalculated tables [Yonsei-Yale isochrones (Kim et al., 
2002) ; ATLAS limb-darkening coefficients (Kurucz, 1992), used when 
T^ff > lOOOOK or logg < 3.5 ; PHOENIX limb-darkening coefficients 
(Claret, 1998, 2000), used when Tg// < lOOOOK and logg > 3.5]. Thus, 
by assuming the masses (Mi^2) of the EB components and the system's 
age, we can look-up the radii {Ri,2), the absolute magnitudes {Magi^2), 
and the limb-darkening coefficients for the binary components. We 
then use these values, as well as the observationally-dctcrmined period 
(P) and combined magnitude out of eclipse {rnagcomh)-, to calculate 
the apparent magnitudes {magi^2) and factional radii (ri,2) of the EB 
components, as follows: 



magi = ma^comfe + 2.5 log l + IQ-MMags-Magi) 

mag2 = magi + {Mag2 — Magi) 



(5) 
(6) 
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a = [G(Mi +M2)(P/2^)^]^/^ ~ 

4.2O6i?0(Mi/M0 + M2/M0)i/3p2/3 (7) 
ri,2 = Ri,2/a (8) 

Besides the epochs of echpses, which can be determined directly from 
the EB hght curve, there are only two additional parameters required 
for us to simulate the light curves of the given pairing: the orbital 
eccentricity (e) and inclination (i). For binaries with short periods 
(< 2 days) and a secondary eclipse precisely half an orbit after the 
primary eclipse, it is reasonable to assume a circular orbit (e = 0). 
Otherwise, one should use the eccentricity derived by an EB model- 
fitting program (we use DEBiL). Finding the inclination robustly is 
more difficult. We employ a bracket search (Press et al., 1992), which 
returns the inclination that produces the best resulting fit. 

To summarize, for every combination of component masses and 
system age of an EB, we can look-up, calculate, or fit all the parame- 
ters needed to simulate its light curve (P, limb-darkening coefficients, 
magi^2) ri,2) epochs of eclipses, e, i), as well as its apparent combined 
color. We systematically iterate through many such combinations. For 
each one we compare the expected light curve with the observations, 
and calculate the reduced chi-squared value {xt)- We also compare each 
observed color {Oc ± ec) with its calculated value (Cc), and combine 
them by defining: score = {wxl+J2c=i[iOc — Cc)/ec]'^)/{w+N). Where 
w is the Xu information weighting. We use w = 1, and assume that 
the smaller the score, the more likely it is that we have chosen the 
correct binary pairing. One can visualized this result using a series of 
score{Mi, M2) contour plots, each with a constant age (e.g. Figure 2.b). 



4. Conclusions 

We have described a novel method for identifying an EB's components 
using only its photometric light curve and combined color. By utiliz- 
ing theoretical isochrones and limb-darkening coefficients, this method 
greatly reduces the EB parameter-space over which one needs to search. 
This approach seeks to estimate the masses, radii and absolute mag- 
nitudes of the components, without spectroscopic data. We described 
two implementations of this method, MECI-express and MECI, which 
enable systematic analyses of datascts consisting of photometric time 
series of large numbers of stars, such as those produced by OGLE, 
MACHO, TrES, HAT, and many others. Such techniques are expected 
to grow in importance with the next generation surveys, such as Pan- 
STARRS (Kaiser et al., 2002) and LSST (Tyson, 2002). 
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Table I. A comparison of the results produced by MECI-express, MECI, and conventional analyses 
with their uncertainties (Lacy et al., 2000, 2002, 2003). The square brackets with numerical values 
indicate the deviation of our results from those of the conventional approach. 





MECI- 


express 




MECI 




Lacy et al. (2000, 2002, 2003) 


Parameter 


Mass 1 


Mass 2 


Mass 1 


Mass 2 


Age 


Mass 1 


Mass 2 


Age 




[M.] 


[M.] 


[M.] 


[XL] 




[M.] 


'M . " 


[Gyi 


FS Mon 


2.9 (AO) 


2.0 (A5) 


1.62 


1.52 


1.4 


1.632 


1.462 


1.6 




[77.7%] 


[36.8%] 


[0.6%] 


[4.1%] 


[0.2] 


±0.012 


±0.010 


±0.3 


WW Cam 


2.0 (A5) 


2.0 (A5) 


1.97 


1.89 


0.5 


1.920 


1.873 


0.5 




[4.2%] 


[6.8%] 


[2.8%] 


[1.0%] 


[0.0] 


±0.013 


±0.018 


±0.1 


BP Vul 


2.0 (A5) 


1.6 (FO) 


1.77 


1.48 


0.7 


1.737 


1.408 


1.0 




[15.1%] 


[13.6%] 


[2.1%] 


[5.4%] 


[0.3] 


±0.015 


±0.009 


±0.2 



We are grateful to Guillermo Torres for many helpful conversations. 
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